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Abstract 

Within the framework of the path-integral approach we study the quantum vor- 
tex creep for the situation where both the Hall and the dissipative dynamics are 
simultaneously present. We calculate the relaxation rate and the crossover tem- 
perature separating the thermal activation and the quantum tunneling processes 
for anisotropic or multilayer superconductors. The results are compared with the 
available experimental data. 
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After the discovery of high temperature superconductors (HTSC) the relax- 
ation of magnetization has been the subject of intensive studies because of the 
possibility of quantum tunneling of vortices in the system. According to the 
Anderson-Kim model, || the vortices can move out of the pinning sites by the ther- 
mal activation process, which proceeds at a rate proportional to exp(— Uo/ksT), 
where Uq is the height of the energy barrier generated by the pinning mechanism. 
This process induces a redistribution of vortices, hence of the current loops related 
with the vortices, causing a change in the magnetization logarithmically with time 

M(t) = Mo[l-Q(T)ln(t/t )], (1) 

where Q(T) = ksT jU§ is the so called magnetic viscosity. Since the rate is 
expected to vanish at zero temperature, the magnetization tends to be unchanged 
as the temperature T approaches zero. However, it has been realized that the 
thermal activation theory fails to explain experimental observations of temperature 
independent magnetic relaxation in HTSC|]lj] at low temperatures. This has led 
to introduce the concept of quantum tunneling of vortices trapped in the pinning 
potential with the collective pinning length L c .@ Since a pinned vortex is in some 
metastable state, the quantum tunneling of the vortex is governed by the frequency 
of the pinning barrier hindering the tunneling process. The rate of the quantum 
tunneling can be thus obtained with the relevant dynamics, whereas the rate of 
the thermal activation process does not depend on particular dynamics. 

Recently, Blatter et al.Jj) studied the quantum vortex tunneling for the case 
where the dissipative term is dominant in the dynamical equation of motion of 
vortex, and many experimental results have been interpreted within this limit. 
On the one hand, Feigel'man et al.0 estimated the low-lying level spacing (~ 
hu>o ~ A 2 ftp) in the vortex core and the transport relaxation time of the charge 
carriers r r at T = by using a linear extrapolation of the normal state resistivity 
p n , and obtained uj T r ^> 1 in HTSC at low temperature which is the superclean 
limit. The criteria distinguishing the dissipative and the Hall vortex motion is the 

2 



magnitude of the quantity WoT n || therefore they proposed that the Hall motion 
should be dominant in clean HTSC Most recently, however, van Dalen et al.p|-|TT| 
observed experimentally that in YBCO or BiSCCO systems, the vortex tunneling 
at low temperatures may occur in an intermediate regime between the purely dis- 
sipative and the superclean Hall tunneling. Prior to the experiments, Stephen |12| 
studied the quantum tunneling of vortex with the dissipative and the Hall terms 
by using the approximate form for the mean-square quantum fluctuation with a 
simple harmonic pinning potential up to a cutoff. As is noted in Ref. since 
the cutoff model potential can be considered to be of the form (q/qo) 2 — (q/qo) n 
with n — > oo where g's are coordinates of the potential, the steep part of the 
model potential underestimates the suppression factor, i.e., the effective action in 
the case of damped systems. So they remarked that a special care should be taken 
in replacing the smooth potential by one which has a discontinuity. 

In this paper, we study the quantum vortex tunneling for the case where the 
dissipative and the Hall dynamics are simultaneously present, considering the pin- 
ning potential barrier generated by disorder with a cubic term|2], |7], [T3J along the 
x axis near criticality, which should better model the real pinning potential. Since 
the tunneling rate takes the form Tq oc exp(—S c \/h), where S c \ is the least ac- 
tion obtained from the classical trajectory, we calculate the classical equations of 
motion in the imaginary time path-integral formalism. When both the dynamics 
are simultaneously considered for the equation of motion, the resultant integral 
equation becomes complicated by the presence of the cubic term in the pinning 
potential, which compels us to rely on numerical methods to obtain the solutions. 

Our considerations will begin by discussing the pancake vortex in the xy plane 
with the length L c along the z axis, where L c is the collective pinning length which 
can be expressed in terms of the mass anisotropy parameter e\ = m/M < 1, 
the coherence length £, and the depairing and critical current densities jo and j c ; 
L c ~ £aCOo/ic) within the weak collective pinning theory. [f2| The supercurrent 
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characterized by the parameter e = (1 —j/jc) 1 ^ 2 is along the y direction. Since the 
length scale in this study is much larger than the size of a vortex core, a vortex 
can be regarded as a point-like object. Neglecting the inertial mass term in the 
equation of motion, Kopnin et al.|| derived the following equation of motion: 

00 

r]u + ail x z = — j x z + F p , (2) 

where 0o — hc/2e is the flux quantum, u the distortion of the vortex from the 
equilibrium position in the xy plane, F p the pinning force, and rj and a are the 
viscous and Hall drag coefficients given by@, || 

(uoT r ) 2 UJ T r 

Here ao = Tchn s , r r = m/ne 2 p n , A = fiVp/TT^, and n s is the density of charge 
carriers, m its effective mass, and vp the Fermi velocity. 

To study the tunneling process at zero temperature, we need to consider the 



imaginary time path integral |T5 1 



J Du x J Du y exp(-S E /h). (4) 

Considering the effect of dissipation on the tunneling, we use the model introduced 
by Caldeira and Leggett |T3| and others 0, [j], [12|] in which dissipation is modeled 
by coupling the vortex degree of freedom to a bath of harmonic oscillators. By 
combining the pinning potential with a cubic term with the Lorentz term in Eq. 
(H) the model potential for the pancake vortex is taken asp], [7|, [TJj 

v p (u x ,uy) — = Y^W l€ ^W ~3 C2 (~R'i' ( ' 

where e = yjl — j/j c , Ci j2 ~ 1, V ~ (4> /4ir\ xy ) 2 , and X xy is the bulk planar 
penetration depth. |16| 



Now, the Euclidean action of this system takes the formfT7|| 

S E = L c /_^r{-^^ % + |[(|) 2 + c ie (|) 2 -?c 2 (|) 3 ] 
1 r°° 

+-/ dnKoir-n^-uin)] 2 }, (6) 



where the so called nonlocal influential function becomes Kq(t) = / °° dujJ(uj) exp(— uj\t\) 
where J{oS) is the spectral density whose form will be discussed later. 

Let us now introduce u x0 , u y0 , and r as the length scales for u x , u y , and r, 
respectively, which are to be set by the following scaling analysis for Eq. 
A comparison of the quadratic potential term with the cubic one in the x direc- 
tion determines the relevant length scale along x\ u x q ~ eRc\/c2- Equating the 
quadratic potential terms in the x and y directions, we obtain the length scale 
for the displacement in the y direction; u y0 ~ e^ 2 Rcl /2 /c 2 . The t ime scale for 
the imaginary time, tq ~ aR 2 /(Voy/cie), is set by comparing the Hall term with 
the quadratic potential in the y direction. Keeping these points in mind, it is 
convenient to use the dimensionless variables 

t cieR, ,ct /2 e 3 / 2 R, ,V2R 2 a , 
u x = (^— K, u y = ( — — )u y , r = (— — — )r. (7) 

The Ohmic dissipation^], |7|, |T2|, Q in Eq. (Q), which is obtained if J{uS) = t]uj 
after elimination of the oscillators, leads to the influential function 

K (t) = r//27r|r| 2 . (8) 

Then, the Euclidean action is simplified as 



du x _ 1 _ 2 1 _ 2 1 



s b = (^-^)(L c a Q R 2 y/ 2 I, (9) 

4c 2 



f°° 7-r • CLU X _ I_ 2 I_ 2 

= J-00 ^ mi lF Uv+ 2 Uy+ 2 Ux ~6 Ux 

+ l m f°° df J^(r) -^(ri)] 2 + ci€[tZ y (f) -%(fi)] 2 (]0) 



where the dimensionless Hall and dissipation coefficient are given by 

1 (^QT r ) 2 V2 U T r 

Ql ~ v^I + Ktv) 2 ' m ~ 27r v ^l + (^ r r ) 2 ' 1 J 

Let us briefly consider the case of a pure Hall type which is relevant in super- 
clean superconductors. In this limit the integral ( [TOD does not include 771 -term. 
Noting that I ~ a\ j^dfu c y du c x / df = a±x (numerical quantity) where u x and 



Uy are classical trajectories in the Hall limit, the approximate least action takes 
the simple form 

S cl ~ (L c aoi? 2 )e 5/2 a! ~ (L c a R 2 )e^ 2 , (12) 

which agrees up to a numerical factor with the result in Table || obtained by the 
explicit instanton solution. In the case where the dissipation term is dominant in 
the vortex dynamics, || ai-term is not included in the integral (10). Noting that 



I ~ 7/i x (numerical quantity) and rji ~ UQT r /\/e in this limit, S c \ is approximately 
S cl ~ (L c a R 2 )e 5 / 2 Vl ~ (L c a # 2 )eW r ), (13) 
which is also consistent with the result in Table [l] up to a numerical factor. fi~9| 



In the intermediate regime, i.e. when the Hall and the dissipative dynamics are 
simultaneously present, the classical trajectories of u x and u y should be determined 
by the Euler-Lagrange equation derived from Eq. flit]): 

• du v u 2 f°° , ,du x , 1 n 

iai-£+v, x --Z-T h I dr x {—^)- = 0, (14) 

ar 2 J-oo ari n — r 

-iai^- + u y - T]xCie ( cZn(^) _ 1 _ = 0, (15) 
ar J-oo ari Ti — r 

where m x (— t) = w x (f) and u y (—f) = —u y (f). As in the dissipative case, the 
equations ([14]) and ([TBI) are nonlocal in time, which necessiates the Fourier trans- 
formation as follows: 

r Oi^) 2 -I . | ,,_ / x 

— 1— r + 1 + 

1 + 7r?7iCie|a;| 
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V87T 

and 



du x u x {u — ui)u x {ux), (16) 



«vH = — — ^ | (17) 
where w^— c<j) = u x {uj) and tt y (— cu) = —u y (uj). Also, we note that Eq. (|i"0"|) is 
reduced to 
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f 00 , ,l r (aiouY , , 

/ <M-h— ■ n + l + 7rr/i a; m^w 

J-oo 2 1 + 7T77iCie|a;| 

1 f 00 

=M a; (a;)u x (a;i)M a; (a; + ^i)}, 
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after Fourier transform and integrating over the u y (u) variable. 

Even though Eq. ( [T6| ) is a one-dimensional problem with respect to u x , its 
analytic solution becomes complex by the presence of the nonlocal term arising 
from the cubic potential (^j). We have therefore numerically solved the above 
integral equations as follows. 

Let us introduce, for a given function T = ^(w), 

U(u) ee f («i^) 2 +1 + ^1!^) (19) 
1 r°° 

K r (uj) = ^= duxFiu-u-fiFiui). (20) 

V87T J-00 

With the function J-(oj) in the above definitions replaced by the Fourier trans- 
formed classical path u x (u>), the equation (|T6| ) corresponds to 

Cr{w)-Kr{w) = Q. (21) 

Let us now define, for a given function J 7 , 

5 r = lim 1 / du \C r {u) - TZ r {u) \ . (22) 

Sl^oo 1 1 Jo 

If a trial function JF has a set of adjustable parameters {pi}, 5?r becomes a function 
of {p^ and if the set {p*} which minimizes 5 can be found, the function JF with 
the parameter set {p*} will be the desired solution of the integral equation. As 
shown in Table [l], a natural choice for the trial function T for general values of ol\ 
and r}\ would be to linearly combine the two extreme solutions such as 

T{bS) = p\U)j sinh(p 2 ^) + Pz exp(p 4 |o;|), (23) 

where pi,P2,P3, an d Pa are parameters to be controlled in the variational procedure. 
We have carried out minimization of 5 of Eq. ( p2|) with the trial function of 



Eq. (23) by the conjugate gradient method. It turns out that the variational 
method is very successful, covering all ranges of Uqt t for arbitrary values of e. In 
principle, a perfect minimization of 5 should make it to be zero, but due to the 



numerical accuracy in carrying out integrations involved in the procedure (or it 
may be possible that the trial function of the form in Eq. (|23|) cannot perfectly 
minimize 5), our numerical calculation typically gave 8/C < 10~ 3 , where L is the 
average value of Cr(uS) in the integrand of Eq. (f^). Nevertheless, the accuracy is 
far good enough for our purpose of evaluation of the action integral in later stages. 
The accuracy of the variational method has been also checked by comparing the 
solutions by this method with those by the direct iterations of differential equations 
([14]) and ( |l"5|) in the range of loqT t where the solutions by the iteration method are 
available. ||21|| The two solutions always match excellently, which confirms that our 
variational method produces true instanton solutions. 

Typical instanton solutions thus obtained are illustrated in Fig. [I]. The clas- 
sical Euclidean action S c \ is given by a numerical integration of Eq. (|T0|) with 
the instanton solutions. We then have the relaxation rate or the magnetic vis- 
cosity Q(T = 0) which equals h/S c \, in other words, Q(0)/Q = 2^/2/1 where 
Q(0) = Q(T = 0) and Qq 1 = e 5 ^ 2 (im s L c ^ 2 ). Here we have set ci 2 to be 1 and 
replaced R by the coherence length £. The results of evaluation of Q(0)/Q versus 
uj T r for e = 1, 0.1, 0.01, and 0.001 are shown in Fig. |2|. We first point out that our 



Q(0)/Qo is different from the result in Ref. |12| , which reflects the fact that the 
infinitely steep potential which underestimates the suppression factor is carefully 
treated in our model potential. p2| As is noted in Table [l], the approximate value 
of e 5 / 2 Q(0) in units of irn s L c ^ 2 is 15/18 in the Hall limit and 2c 1 / 2 / '(iTujoT r ) in the 
dissipative limit, which agree with the limiting values of the figure. Even though 
Q(0)/Qo(— fc/£>ci) looks going toward infinity as u r r — > 0, it actually does not 
diverge in that limit. In the pure dissipative case, the classical action S c \ vanishes 
as uoqt t — > in our consideration, which is expected from the limiting behavior of 
Eq. (p~3|) . However, in this situation the inertia term not included in this work 
influences the results. By including the massive term, in this limit the approximate 
form of the classical action becomes S c \/h ~ L c \Jm v V^e'l' 1 which is independent 



8 



of u>oT r where m v is the inertia mass of a vortex. 0, Q In general the massive 
term is small in the region except for the limit mentioned before and usually can 
be neglected as compared to either the dissipative or the Hall term. In the inter- 
mediate tunneling regime, our curves display an interesting feature which is the 
clear existence of minimum of Q(0)/Qo for each e. Even at e = 1, there exists a 
shallow yet clear minimum near ujqt t ~ 3.5 (see the inset). As e becomes smaller, 
i.e., as j — > j c , the minimum becomes much more pronounced with its location 
moving toward uj r r « 1 at the same time. The existence of the minimum of Q 
can be understood by noting the following qualitative features of the relaxation 
rate based on the dimensional estimation. || Noting that for the pure Hall motion 
the classical trajectory is dominant in the frequency scale \u\ ~ 1/oti in Eq. ([16]) 
which is also understood by the solution u x (uj) in Table [l|, we find that from Eqs. 
(P) and (|T8|) the correction to the pure Hall classical action Sq by the small dis- 
sipation rji is of the order rji/ai, that is, the corrected classical action becomes 
(l + rji/ai)SQ ~ (1 + \/ujQT r )SQ . The relaxation rate is then Q ~ uor r /(l +u T r ), 
which increases monotonically as a function of ooqt t in the large ujqT t regime. The 
fact that the classical action increases (or the corresponding relaxation rate de- 
creases) by inclusion of the dissipation is well known in the macroscopic quantum 
tunneling problems, which is also physically clear in the sense that the dominant 
(Hall) motion is hindered by the frictional force, resulting in decreasing the relax- 
ation rate. In the opposite limit, on the other hand, the correction to the pure 
dissipative action Sq by the small Hall contribution a± is of the order (cti/771) 2 . 
The corrected classical action is then (1 + (ax/r^) 2 )^ ~ uj r r (l + (u T r ) 2 ), leading 
to the relaxation rate Q ~ l/[u)oT r (l + (u} T r ) 2 )]. Q now decreases monotonically 
as a function of uoT r in the small u>oT r regime. Therefore, it is expected that there 
exists a minimum of Q in the intermediate range of Uqt t because the relaxation 
rate decreases monotonically for small ujqt t but increases monotonically for large 
uoqt t . This suggests us to predict observation of minimum in an experiment. 
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As is shown in Fig. ^, our results are valid in any value of e < 1. It is well 
known that the tunneling rates get larger for smaller e since the height and the 
width of the barrier are proportional to the parameter e. However, there might 
exist the quantum tunneling even in the absence of the external current, i.e., e = 
1. Actually the existence of the quantum tunneling is determined by estimating 
the magnitude of the quantum tunneling rate Tq. If we apply our result of e = 
1 to the existing experimental data, [16] where Q(0)/Qo ~ 2.3 and Qq 1 ~ 78 



in the YBa2Cu30 7 /PrBa 2 Cu307 multilayer systems|| and Q(0)/Qo ~ 2.0 and 
Qq 1 ~ 66 in the BiSr 2 CaCu20s single crystal, flToH we obtain S c \/h ~ 33 for both 



systems. This gives Tq « c<jexp(— S c \/h) « 4.66 x 10 3 (4.66 x 10 4 )sec 1 for the 
attempt frequency w(w V e l l 2 / '(£ 2 n s )) of order 10 12 (lO^sec" 1 for YBCO/PBCO 
(BiSCCO), which is a resonable value to ensure the observable rate of tunneling. 
Also, using the value of Q(0)/Q for both systerms, we get u r r « 0.29, which 
corresponds to the Hall angle 0# = arctan(ai/?yi) = arctan(a>oT r ) ~ 16° for the 
former, and Uqt t ~ 0.37 (0# ~ 20°) for the latter. This implies that for those 
systems the vortex tunneling at low temperatures occurs in an intermediate regime 
between the purely dissipative tunneling and the superclean Hall tunneling. It 
should be also noted that the quantum creep in the oxygen deficient YBa2CusO x 
films [0] occurs in the intermediate regime as well, where the Hall angle depends 
on the oxygen content. 

In the thermal activation regime, the classical solutions u x (f) and u y (f) do not 
depend on rp4| and the integration range in Eq. flTOf) reduces to between —(3/2 



and (5/2 where (3 = (3U(^[c~^V§/ \Z2R 2 a ). Then we obtain the classical action 
S c \ = St = ftfiUo, where Uq = L c Voe 3 /6, and the corresponding escape rate 

T T cx exp(-S T /h) = eM-U Q /k B T), (24) 

which is the Boltzmann formula representing a pure thermal activation. Com- 
paring T T with rg(cx exp(— S c \/h)), we have the crossover temperature k B T ~ 
Ae 1 l 2 Q(Q) / Qq from the thermal to quantum regime where A = 0Q/(967r 3 A 2 £ 2 n s ) 
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and find that the quantum process dominates at T < Tq. Using £ = 1.5nm (2.0nm), 
A = 150nm (250nm),]25| and n s = 5 x 10 27 /m 3 (3.5 x 10 27 /m 3 ) for YBCO mul- 
tilayer systems (BiSCCO single crystal), for e = 1 we get T ~ 9.48 K (2.36 K) 
which agrees reasonably with the experimental results. 

In conclusion, we have studied the quantum vortex creep for the case where 
the Hall and the dissipative dynamics are simultaneously present by using the 
standard instanton method. We have introduced the cubic potential term to better 
model the pinning potential and solved the resultant integral equations numerically. 
In a comparison with available experimental data, our results indicate that the 
quantum vortex creep may occur in an intermediate regime between the Hall and 
the dissipative regimes for highly anisotropic and multilayer superconductors. It 
has been also found that there exists a minimum in a scan of the relaxation rate by 
u T r for every e we have studied, whose observation we predict in an experiment 
with nonzero external currents. The temperature is set to zero in this study, but 
an extension to finite temperatures will be discussed elsewhere. It will be also 
interesting to apply our approach to systems with columnar defects produced by 
irradiation with heavy ions. 

This work was supported by Research-Fund, Korea Research Foundation 1996 
and by the Ministry of Information and Communications of Korea. 
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Table 1: Instantons in the r space and in the u space, and the corresponding 
actions in the Hall limit (u> T r ^> 1) and the dissipative limit s(co> TV <C 1). Here 
ch = 18c^ 2 /15c2 and cd = -kc\/2c\. 
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Figure 1: Typical instanton solutions: u x (r) (top) and —iu y {j) (bottom) for e = 1 
when u> T r 3> 1 (a), u r r = 1 (b), and u r r = 0.35 (c). 



Figure 2: Dependence of the quantum vortex creep relaxation rate Q{0)/Qo on 
u Q T r at zero temperature, for e = 0.001 (a), 0.01 (b), 0.1 (c), and 1 (d). Also is 
shown, to logarithmic accuracy, the relaxation curve (e) from Ref. [12[ (Eqs. (16) 
and (20)). Note that the approximate value of Q(0)/Qo becomes 15/18 in the Hall 
limit (ujQT r ^> 1) and 2e 1/ ' 2 /( 7ru; o r r) in the dissipative limit (co T r 1). In the inset 
the plot (a) for e = 1 is magnified. 
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